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The use of single molecules to design electronic devices is an extremely chal- 
lenging and fundamentally different approach to further downsizing electronic 
circuits. Two-terminal molecular devices such as diodes were first predicted |T] 
and, more recently, measured experimentally |2]. The addition of a gate then 
enabled the study of molecular transistors [3H5]. In general terms, in order to 
increase data processing capabilities, one may not only consider the electron's 
charge but also its spin |6l This concept has been pioneered in giant mag- 
netoresistance (GMR) junctions that consist of thin metallic films [S], Pj. Spin 
transport across molecules, i. e. Molecular Spintronics remains, however, a chal- 
lenging endeavor. As an important first step in this field, we have performed an 
experimental and theoretical study on spin transport across a molecular GMR 
junction consisting of two ferromagnetic electrodes bridged by a single hydro- 
gen phthalocyanine (H2PC) molecule. We observe that even though H2PC in 
itself is nonmagnetic, incorporating it into a molecular junction can enhance the 
magnetoresistance by one order of magnitude to 52%. 

Scanning Tunneling Microscopy (STM) has proven to be a powerful and versatile tool for 
studying electron transport properties of single molecules, be it by continually addressing 
individual molecules on a one by one level [T0Hl3], or in an automated break junction 
protocol [111115]. In this work, we introduce spin-polarized STM (Sp-STM) to measure the 
magnetoresistance of single molecules, employing spin-polarized electrodes. Experiments 
were carried out in a home-built STM working in ultra high vacuum at 4 K [TB] on individual 
H2PC molecules (C32H18N8) sandwiched between Co-coated W tips and ferromagnetic Co 
nano-islands on Cu(lll) single crystals. 

In Fig. [1] we present the STM topography of the sample after depositing H2PC molecules, 
that identify themselves by their four aromatic isoindole (BzPy) side groups [131 HI] • We 
note that the Co nano-islands exhibit a spontaneous out-of-plane magnetization due to a 
strong surface anisotropy [H]. By using Co-coated tips (10 monolayers) with an out-of-plane 
magnetization, we can use the Sp-STM technique's sensitivity to the spin-polarized density 
of states [19] to determine the relative orientation (parallel:P or antiparallel:AP) of the 
magnetization of individual islands relative to that of the tip. Fig. [T)d shows the differential 
conductance (d//dV^) curves measured atop two islands with P and AP alignments with 
respect to the tip. Particularly large differences in the spectra are found at -350 meV, which 
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corresponds to the surface state of Co [20]. This difference allows us to detect the local 
magnetization direction of Co islands by recording maps of the local differential conductance 
at this bias voltage (see Fig. [l^). The differential tunneling magnetoresistance (differential 
TMR), defined as the difference in conductance divided by the smaller conductance of the 
tunneling junction formed by the tip and sample, is strongly energy-dependent as depicted 
in Fig. We note here that the differential TMR for this Co/vacuum/Co junction is only 
~ 5 % at low bias voltage. 

To contact single molecules, we positioned the STM-tip above the aromatic side groups 
of a H2PC molecule, opened the feedback loop and decreased the tip-to-molecule distance 
(^ lA/s) while measuring the tunneling current. We thus obtain conductance-distance 
curves [T014T2] . Since the observed conductances of phthalocyanine molecules are rather 
high [13], a small voltage (10 mV) has to be used to avoid the thermal disintegration of 
the molecules. We observe an exponential increase of the conductance that reflects the 
decreasing tunneling barrier width (see Fig. |2^). Below a certain tip-to-surface separation 
— typically 3-4 A — the conductance abruptly increases, which indicates a sudden change 
of the junction geometry. Similar jumps have been seen, e. g., for alkanedithiol molecular 
wires [TT]. We have recently shown that, for phthalocyanine molecules, this corresponds to 
a lifting of the flat H2PC molecule's aromatic group as it contacts the tip fT3]- 

The conductance after the molecular jump-to-contact depends only slightly on the dis- 
tance and encompasses both the transport across the molecule Gmoi and direct tunneling 
between the tip and the sample Gtun [Hj (see Fig. 2c). This conductance G in molecular 
contact mode depends markedly on whether the underlying island is of P or AP type. In 
the particular measurement of Fig. [2^, we find Gp = 0.30 6*0 and Gap = 0.17 Go) measured 
at 10 mV bias, where Gq = ^ is the quantum of conductance. We may eliminate a sizeable 
proportion of the spin-dependent direct (tip-to-island) tunneling contribution to the con- 
ductance G by subtracting the measured conductance Gtun before the jump from that after 

Gcont) Gmol ~ Gcont Gtun- 

To ensure identical tip conditions when quantitatively comparing P and AP conductances, 
measurements on two Co islands of P and AP type were performed in the same scan. This 
approach further eliminates magnetostriction of the two ferromagnetic electrodes. Each 
such measurement was in turn repeated several hundred times and the distribution of the 
P and AP conductances is depicted in Fig. [2ji. The width of the conductance distribution 
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- which underscores the noise of the conductance measurement and a variation in contact 
geometries - is relatively small when compared to previous work [TT]; the histogram of the 
measurements on nearly 800 P and AP junctions clearly reveals a difference in conductance 
between the P and AP junctions. Gaussian fits were used to determine the GMR from the 
measurement statistics. We find Gp = (0.242 ± 0.007) Go and Gap = (0.160 ± 0.005) Gq. 
These values result in an optimistic GMR ratio of 

GMR= ^^~^-^^ = (51±9)%. 
Gap 

Surprisingly, the GMR ratio obtained at V = 10 mV is one order of magnitude larger than 
the differential TMR found for direct tunneling between the tip and the Co surface. 

To understand what causes this large value of GMR, we have performed transport cal- 
culations based on density functional theory (DFT) employing the nonequilibrium Green's 
function (NEGF) formalism and the TURBOMOLE package [21]; details may be found in 
the Supplementary Information. 

To find the atomic structure of the molecular junction, we first performed a geometry 
optimization for H2PC on the Co(lll) surface, represented by a cluster with 65 atoms. Our 
analysis suggests that H2PC adsorbs preferentially in the bridge position onto Co(lll) (see 
Fig. [3^), owing to a binding energy —8.17 eV that is larger than that found in either the 
hollow site position (—8.06 eV) or the atop site position (—7.45 eV), consistent with earlier 
findings [22l|23]. 

We have calculated spin-polarized transport in the linear response at low bias voltage 
for the two junction geometries schematized in Fig. |3^,b) before ("fiat") and after ("con- 
tact") the jump-to-contact. To establish the latter configuration, a free H2PC molecule has 
been bent along the low energy vibrational eigenmode with frequency u = 6.4 meV [T3] . 
Within the NEGF calculations, the magnetization direction (P or AP relative to the tip) 
of each Co-cluster is a control parameter. We have ascertained that, in our study, the elec- 
tronic structure of the Co surface is properly reproduced with, in particular, an exchange 
splitting of the li-states of 1.8 eV (see Supplementary Information). This leads to a mag- 
netic moment m(Co) ~ 1.65 pei' surface Co atom, in agreement with previously reported 
calculations 

On a qualitative level, our transport calculations (see Fig. ^) reproduce very well our 
experimental findings (see Fig. [2^). We confirm the exponentially increasing conductance in 
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the tunneling regime, Gtunid) ~ e"'''^, for which the distance d between the two electrodes 
is still large. The slope is independent of the relative alignment of electrode magnetizations 
and the computational value = 1.87 A^^ (i. e. the work function = 3.24 eV) is 

in agreement with experiment (/3 = 1.9 ± 0.3 A; W = 3.2 eV). 

Once contact between the tip and the molecule has been established, G'cont('^) varies much 
more weakly with the contact distance d just as observed in the experiment. It is, however, 
still sensitive to the relative orientation of the magnetization of the electrodes. We find that 
GAp{d) is always much lower than Gp{d). For a quantitative comparison with experiment, 
we consider the GMR*^''° ratio at the distance d for which the ratio r = Gtun/ Ccont matches 
the value r ^ 4 found experimentally. We thus find that GMRtheo ~ 65 % and is only weakly 
dependent on d. 

We now discuss the conduction mechanism from the substrate onto the molecule and 
across the molecular contact established between the tip and the aromatic group, and its 
impact on the large GMR measured. Pc molecules are characterized by an energetically 
isolated highest occupied molecular orbital (HOMO) and a nearly doubly degenerate lowest 
unoccupied molecular orbital (LUMO) [25]. We note that the HOMO* levels corresponding 
to the aromatic group hybridize only very weakly, with almost no amplitude on the bridging 
nitrogen (Nb)4. By contrast, the LUMO states are located on two out of the four aromatic 
groups, with a strong hybridization to all N atoms forming the inner macro-cycle (see Sup- 
plementary Information). Since the N bond to Co includes states at the Fermi energy Ep 
[13], transport should occur via the quasi-degenerate LUMO level. We confirm this fact by 
examining in Fig. [3]i) the transmission probability per spin direction T'^^i(E) at the Fermi 
energy across a junction of P type. We find that Gp near Ep is indeed weighted by a peak 
centered slightly above Ep that underscores transmission through the LUMO level. As we 
discuss in the Supplementary Information, the larger density of Co minority states at Ep 
results, through the N-Co bond, in a larger efficiency of LUMO hybridization, and thus of 
LUMO broadening, in the spin J, channel. 

This difference in LUMO broadening for the two spin channels has a direct impact 
on the GMR measured across the molecular junction. Indeed, the conductance across 
a single level, here the LUMO, generically takes on the Breit-Wigner form [26j G ~ 
psubstrateptip/(^^^^^j^^ - Epf + (F^^'^'^^rate ^ r*'P)V4), which couslders the energy separation 

between the LUMO and Ep , as well as the LUMO broadenings (inverse lifetimes) P'^"^'^*'''^*'^ 
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and r*'P due to hybridization to the substrate and tip, respectively. Each is in turn spht 
into r™™^™'^^) depending on the spin channel considered. Because transport is off-resonant, 
i.e. I^LUMO - ^fI > r-^'^'^^j, we have Gp ^ (r'^^^r'"'^ + r'^^jr'^'^j)/(ELUMO - ^f)^ while 
Gap ~ 2r"^''^r"^^V(^LUMO - ^f)^- Introducing the ratio g = T^^^i /V^^'^ we thus find 

GMR«'r--»-r-.-y^a-ri! 

2pmmpmaj 2^ 

This simple formula implies two important rules of thumb for spin-polarized transport off- 
resonance across a molecule. First, the GMR is insensitive to the precise location of the 
resonance energy. Second, it is mainly indicative of the ratio g of minority and majority 
molecular orbital broadenings due to hybridization. Relative to the small differential TMR 
ratio found between the Co substrate and tip, the larger GMR ratio measured experimentally 
can be explained theoretically by the above ratio g. The imbalance within the two spin 
channels of the hybridization-induced level broadening to the molecular orbital responsible 
for transport promotes a high GMR across the molecule. 

In conclusion, we have experimentally demonstrated a GMR of over 50 % in single 
molecule junctions that is much larger than the value of differential TMR found without the 
presence of molecules in the junction. Such a high value is caused by a strong hybridization 
of the molecular LUMO, responsible for transport, with minority states of the two metallic 
electrodes. Such a selective hybridization thus leads to a spin filtering effect that could be 
generic to all molecular junctions. 

We express our gratitude to O. Hampe, J. Kortus, K. Fink, S. Boukari, Xi Chen, M. 
Alouani, R. Mattana, J. van Ruitenbeek and P. Seneor for useful communications and ac- 
knowledge support by the DFG (WU 349/3-1 and SPP1243), the Alexander von Humboldt 
foundation and the ANR (ANR-06-NANO-033-01). 

Methods 

The Cu(lll) crystal was cleaned thanks to several cycles of Ar+ sputtering and annealing. 
The molecules were evaporated in situ from a Knutsen cell heated to ~ 500 K. During the 
deposition process we keep the sample at 270 K to reduce thermal diffusion of the deposited 
molecules. d//dV^ curves were measured on the bare islands with the lock-in technique. 

DFT-based transport calculations were carried out with a homemade code building upon 
the NEGF formalism and the TURBOMOLE package |21]. Our implementation enables us 
to perform transport simulations with free boundary conditions, which for the present case 
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have been extended to account for the spin-polarized electronic structure of the magnetic 
electrodes (for further details, see Supplementary Information). The gradient corrected ap- 
proximation (GGA) DFT-energy has been amended by empirical corrections [27j to account 
for dispersive van der Waals interactions between the molecule and the surface. 
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FIG. 1: a) Topographic image of H2PC molecules adsorbed onto two Co islands on the Cu(lll) 
surface. The colour code displays the measured dl/dV map at -310 mV. One can distinguish 
between the two islands species with a magnetization parallel (P; in yellow) and antiparallel (AP; 
in red) with respect to the tip magnetization, b) dl/dV spectra taken on two islands of P and AP 
type, which clearly reveal spin-polarized states below the Fermi edge, c) Optimistic TMR ratio 
calculated from the dl/dV spectra. The highest value is measured at around -350 meV, and is 
used to distinguish between the two islands species. 
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FIG. 2: a) A typical set of conductance-distance curves measured atop a H2PC molecule adsorbed 
onto P and AP magnetized islands with a constant tunneling voltage of 10 mV (Go = As the 
tip approaches the molecule, the tunnel barrier width decreases, so that the conductance increases 
exponentially (see panel b)). Below a certain tip-to-surface separation — typically 3-4 A — the 
conductance abruptly increases as the molecule jumps into contact (see panel c), and then varies 
only slightly upon further reducing the distance. Transport across the contacted molecule reflects 
both direct tunneling between the tip and the surface and conduction across the molecule (see 
panel c)). d) Histogram of the corrected molecular conductances (390 times parallel / 384 times 
antiparallel) . A Gaussian fit is used to determine the statistical conductance in the P and AP 
configurations, and thus the GMR ratio. 
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SUPPLEMENTARY INFORMATION: 



Magnetoresistance through a single molecule 

Stefan Schmaus, Alexei Bagrets, Yasmine Nahas, Toyo K. Yamada 

Annika Bork, Martin Bowen, Eric Beaurepaire, Ferdinand Evers, and Wulf Wulfhekel 

1. Geometry optimization 

We used the quantum chemistry package TURBOMOLE [Ref.l] to find the atomic struc- 
ture of molecular junctions. We have performed standard geometry optimization for H2PC 
on Cu(lll) surface, represented by a cluster with 65 atoms. (The use of Co(lll)-clusters for 
structure optimization is impaired by the very large number of energetically almost degen- 
erate spin multiplets.) The gradient corrected approximation (GGA [Ref.2]) DFT-energy 
has been amended by empirical corrections [Ref.3] to account for dispersive van der Waals 
interactions between the molecule and the surface. Our analysis suggests, that H2PC prefers 
the "bridge" position (see the paper. Figure 3a) with binding energy —8.17 eV, against the 
"hollow" site position (binding energy —8.06 eV) and "atop" site position (binding energy 
—7.45 eV), which is consistent with earlier findings [Refs.4,5]. 

2. Electronic structure and transport calculations through H2PC 

Density functional theory (DFT) based electronic structure and transport calculations 
through H2PC molecular junctions have been performed within the non-equilibrium Green's 
function (NEGF) approach as implemented in a homemade simulation code [Ref.6] interfaced 
to the quantum chemistry package TURBOMOLE [Ref.l]. The atomic configuration of the 
"extended molecule" used to simulate a bottle-neck of the molecular junction is shown in 
Suppl. Fig. [1^ ("contact" regime): H2PC is bound to the two Co(lll) clusters with 51 and 
19 atoms, representing the Co surface and the Co STM-tip, respectively. To establish this 
configuration a free H2PC molecule has been bent along the low energy vibrational eigenmode 
with frequency u = 6.4 meV. A similar atomic configuration, with H2PC bound to the 
surface only, has been used for transport simulations in the tunneling regime (see the paper. 
Figure 3b). The generalized gradient approximation (GGA, BP86 functional [Ref.2]) and a 
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contracted Gaussian-type split-valence basis set with polarization functions (SVP) [Ref.7] 
have been employed for calculations. 

First, a closed-shell (nonmagnetic) solution for the "extended molecule" comprising 2156 
electrons has been found. To account for infinite reservoirs with spin-polarized electrons, the 
spin-dependent (a = ±1/2) local self-energies, ^^"'^'^^'^'^''^'''(x, x') = [A + crAex — i'j] 5(x — x'), 
have been ascribed to the outermost boundaries of the simulation cluster (dark gray atoms 
at Suppl. Fig. [l^). Here, a parameter Acx = 1.6 eV accounts for exchange splitting of the 
bulk Co (i-states [Ref.8]. A freedom to choose a sign of Aex independently for the "surface" 
and " STM-tip" clusters allows for two solutions: with parallel and antiparallel alignment of 
electrodes' magnetizations. The non-equilibrium Green's function formalism is employed to 
evaluate the charge- and spin-density matrices in the presence of open boundaries insuring 
a charge neutrality within the "extended molecule". The density matrices are given back 
to TURBOMOLE to find a modified set of Kohn-Sham orbitals, with a cycle to be repeated 
unless the self-consistent solution is reached. For the given value of the level broadening, 7 
( = 1.36 eV in present calculations), the contribution A to the real piece of the self-energy 
has been defined by imposing the condition of spurious charge accumulation to be absent at 
the cluster's boundaries. Further details of the computational approach will be published 
elsewhere [Ref.9]. 

As an example, details of electronic structure for the molecular junction in the case of 
parallel alignment of magnetizations are presented in Suppl. Fig. [T]3,d. The electronic states 
at the Co-surface are spin-polarized with exchange splitting ^ 1.8 eV (Suppl. Fig.[l|i), giving 
rise to an average magnetic moment ^ 1.65fiB per a surface Co atom. Inspecting the local 
density of states (DOS) at H2PC, Suppl. Fig. we observe the majority spin resonance 
above the Fermi level (Ep) to be identified with the quasi-degenerate LUMO, which has a 
significant weight on the nitrogen atoms (cf. Suppl. Fig. [i|d). In contrast, no density is seen 
on nitrogens for the majority spin HOMO resonance positioned ~ 1.2 eV below Ep that 
reflects the stucture of the HOMO (Suppl. Fig. [i]d). Furthermore, the minority spin H2PC 
LUMO level evolves into a broad peak due to a hybridization with the minority spin Co 
states, which density near Ep significantly exceeds the ones for majority spin electrons. The 
asymmetry in the LUMO level broadenings, F™"^ > r™^j, gives rise to a magnetoresistance 
effect (see text of the paper for further details). 
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Suppl. Fig. 1: (a) Atomic cluster, representing a bottle-neck of the H2PC molecular junction, used 
for the electronic structure and transport calculations; outermost boundaries of the Co clusters 
(dark gray) are subject to the absorbing boundary conditions modeled by the self-energy Eg-, 
(b) Frontier molecular orbitals of H2PC (D2/1 symmetry): an A„ HOMO and a quasi-degenerate 
LUMO doublet, B23 and 63^. (c) and (d): spin-polarized local density of states at H2PC and 
Co(lll) surface, respectively. 
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